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Abstract 

Using the exact propagators in a constant magnetic field, the neutrino self-energy has been 
calculated to all orders in the field strength B within the minimal extension of the Weinberg- 
Salam model with massive Dirac neutrinos. A simple and very accurate formula for the self- 
energy is obtained, that is valid for < i? ^ m^/e and for neutrino transverse momentum 
to the magnetic field p± <^ mw- I discuss the implications of this finding to the dispersion 
of massless neutrinos in vacuum and in a charge-symmetric medium, and to the magnetic field 
induced resonance transitions of massive neutrinos inside supernovae and magnetars, and calculate 
the neutrino magnetic moment. 
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I. INTRODUCTION 



The study of creation, propagation, energy loss, and absorption of neutrinos in magnetic 
field is important in many astrophysical contexts and in early cosmology The neutrino 
self-energy and dispersion relation are modified in magnetized media and in a vacuum with 
magnetic field, and such modifications have been studied extensively in the literature 2], 
y, y, y, y, Q, y, yi. There is & natural scale for the magnetic field strength required to 
significantly impact quantum processes, and it can be expressed in terms of the electron 
mass me and the elementary charge e as = ml/e ~ 4.41 x 10^^ G. Large magnetic 
fields are present in a variety of astrophysical sites like supernovae, neutron stars and white 
dwarfs, and fields as large as B^. or larger can arise in supernovae explosions or coalescing 
neutron stars. The remnants of such astrophysical cataclysms are magnetars, young neutron 
stars with magnetic fields 10^^ — 10^^ G [lO|, lul, 112] • It has been suggested that during the 



electroweak phase transition local magnetic fields much stronger than those of a magnetar 
could have existed, with field strength as high as 10^^ — 10^^ G jisl, Q, Q. Situations 
where even stronger magnetic fields could exist in extreme astrophysical and cosmological 
environments are possible. While neutrinos might rarely encounter magnetic fields larger 
than 10^^ G, many situations arise where an abundant production of neutrinos occurs in 
astrophysical sites such as supernovae, neutron stars, white dwarfs and magnetars where 
magnetic field strengths can be at or around B„. A literature search reveals that, while 
calculations of the neutrino self-energy in magnetic field have a long history 
only in the last few years has this matter been partially settled with the paper by Kuznetsov 
et al. 6| who, using an expansion for the ly-propagator where only the lowest order terms 
are retained, find the correct asymptotic values of the self-energy for neutrino transverse 
momentum pj_ ^ niw in the case of weak magnetic field eB <^ and of moderate 
magnetic field ml <^ eB ^ m'^. We still do not know the self-energy or dispersion relation 
of neutrinos with low transverse momentum for magnetic field strengths that are not much 
smaller and not much bigger than B^, and therefore a calculation of the neutrino self-energy 
and dispersion relation for pj_ <^ m\y that is valid for magnetic field strengths covering the 
whole range < eB <^ is needed. 



In this paper 1 use Schwinger's proper time method 



16| to calculate the neutrino self- 



energy in homogeneous magnetic fields within the minimally extended Standard Model 
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of electroweak interactions with massive Dirac neutrinos. Using the exact PF-propagator 
I obtain a simple and very accurate analytic form of the self-energy that is valid for 
< eB <^ and when the neutrino transverse momentum to the magnetic field is 

p± ^ iTiw, and I show that, in the case of weak and moderate fields, my result agrees 
with the results of Ref. |6i]. I discuss the implications of this finding to the dispersion of 
massless and massive neutrinos in a plasma and in a vacuum with magnetic field and use it 
to calculate the neutrino magnetic moment. 

In Section [TT]the notation for the fermion, gauge boson and scalar propagators in magnetic 



field ISll is reviewed and the one-loop neutrino self-energy is set up in the framework of the 

nn 

minimal extension of the Standard Model [g, [ITJ. In Section HIT] I calculate the self-energy 
and obtain simple analytic expressions for all its terms. An extended discussion of the 
implications of my results and the conclusions are in Section IIVI 

II. PROPAGATORS AND NEUTRINO SELF-ENERGY IN A CONSTANT MAG- 
NETIC FIELD 

The metric used in this paper is g'^" = diag(+l, — 1, — 1, — 1) and the 2-axis points in 
the direction of the constant magnetic field B. Therefore the electromagnetic field strength 
tensor F^'^ has only two non-vanishing components F"^^ = —F^^ = B. 

For the purpose of this work, it would seem convenient to work in the unitary gauge 
where the unphysical scalars disappear. However, the VT-propagator is quite cumbersome 
in this gauge, and I prefer to work in the Feynman gauge, where the W^-propagator has a 



much simpler expression. The following expressions for the charged lepton S{x', x") 16l. Il9|. 

n 

ly-boson G^^ {x\x") and scalar propagators D[x\x") \b\ in a constant magnetic field have 
been written using Schwinger's proper time method: 

x") = n{x', x") J ^e-^'=-(^'-^")^(A;) , (1) 
= n{x',x") I ^e-^^-(^'-^")G^^(fc) , (2) 

D{x\ x") = nix', x") I ^e-^'=-(^'-^")D(fc) , (3) 

and, in the Feynman gauge, the translationally invariant parts of the propagators are 
ds 



S{k) = I ^exp 

JO cos eBs 



-is I m-, — Km — A; I 

" eBs 



cos ei?s 



, (4) 
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Jo 
D{k) 



°° ds 



exp 



cose-Bs 

ds 



tan eBs 



eBs 



COS eBs 



exp 



—IS I m^r — k,, — k' 



tan C-Bs 
- eBs 



(5) 
(6) 



where — e and are the charge and mass of the charged lepton i, and mw is the VT-mass. 
It is convenient to use the notation 



off = (a°,0,0,a^), a^ = (0,a^a^O) 



(7) 



and 



(a6)|| =a°6°-a^^ (a6)^ = -a^ 6^ - 6^ (8) 
for arbitrary four-vectors a and b. Using this notation I write the metric tensor as 



with 



9"" = 9^ + 



(9) 



(10) 



where (p is the dimensionless electromagnetic field tensor normalized to B and (p is its dual 

(11) 



The 4x4 matrix that appears in the charged lepton propagator (jll), can be written in 
terms of if as 



= 2[V,7l = -i^invi) 



(12) 



where the Lorentz indices of vectors and tensors within parentheses are contracted, e.g. 
(7^97) = •y^ip^'^'y^. When writing the ly-propagator ([5]), I use the notation 



(^^2eFs\^ ""^ = -^f COS i2eBs) + if"'' sin {2eBs) . 



(13) 



I choose the electromagnetic vector potential to be An = —^^F^yX" and therefore the phase 
factor which appears in Eqs. ([T]), ([2]), 



is given by pLS] 



n{x\x") = exp ( -ilx^F'^V; 



(14) 



The neutrino self-energy operator S(j9) is defined in terms of the invariant amplitude for 
the transition Ui — > ui 

M{yi ^ Vi) = -P(p)S(p)z/(p). (15) 
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Perturbatively, the self-energy operator in the Feynman gauge corresponds to the sum of 
two diagrams, a bubble diagram with the gauge boson and a bubble diagram with the scalar 



E{p) = Ew{p) + ^M- (16) 

The translationally non- invariant phase factors Q{x', x") are identical for all propagators 
and the product of phase factors in the two-vertex loop is 

Q{x',x")Q{x",x') = 1 (17) 

therefore, within the minimally extended version of the standard model of electroweak in- 
teractions with an SU (2)-singlet right-handed neutrino, the two bubble diagrams can be 
written as 

m 

^w{p) = -^jRlc./ -^,S{p-k)G^^{k)^pL, (18) 
S*(p) = ^^^:^["^^^ ~ "^'^^l J '(2ttY^^^ " k)D{k)[mgL - m^R], (19) 



where g is the SU{2) coupling constant, L = |(1 — 75) and R = |(1 + 75) are the left- 
handed and right-handed projectors and neutrino mixing is allowed by taking a nondiagonal 
neutrino mass matrix m^, in Eq. (fT9|) . 

III. CALCULATION OF THE SELF-ENERGY 

Inserting the expression for the propagators from Eqs. (jl]) and ([5]) into the self-energy, I 
write Svi/(p) as 

2 J (2tc) J cos 2:1 J COSZ2 



Ria 

where 



COS^i 



^g,a _ ^e'eFs^f^-^^^L (20) 



q = p — k , zi = eBsi , Z2 = eBs2 ■ (21) 



I do the straightforward 7-algebra, change variables from Sj to Zi, translate the k variables 
of integration as follows 

(^il , k^) - ih + -^P\\ , ki. + p^) (22) 

Z1 + Z2 tan 2:1 -I- tan 2:2 
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and, finally, perform the four gaussian integrals over the shifted variables k. The result is: 

-■2 ""^ ""^ dzidz2 
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(47r)2 7 0^0 {zi + Z2) sin(zi + Z2) 

Z2 



Z1+Z2' 



ia-i{zx+2z2) _^ Sin2;2 



sin(2i ^ Z2)' 



L + (c.t. 



where 



V 



Z\Z2 2 sin2;isin22 2 



(2:1 + 2;2j" " sm \^zx + -22 
and the appropriate counter-terms (c.t.) are defined such that 



(c.t.) 



w 



- i> 






B=0,j/=0 ^ 




B=0,j/=0 



(23) 



(24) 



(25) 



Next it is convenient to change integration variables from (zi, Z2) to (t, m) defined by 



u 



and 



Z2 



TU 



(26) 



and to perform a clockwise rotation in the complex plane so that r = —iz. The result is 

.2 i-i 



Sh/(p) 



(47r)2 7 sinh zJ 



R 



u cosh(z + zuj-j^w + u sinh(2; + zu)a-i'^\\ + 



sinh zu 



where I use ""^ 



pjl ^ —p\_ and introduce the parameters r] 
i = p^/m"^ and the function 

sinh 2;m sinh (2; — zu 



sinh z 

eB/m 



2 

1^5 



w(-u, z) = u{l — u) 



L + {cX.)w (27) 
\i = mj/m^^ and 

(28) 



z sinh z 

At this point we must analyze the role played by the neutrino transverse momentum. A 
detailed analysis of the role of p5_ is done in Ref. p^] where the neutrino self-energy is 
calculated in terms of integrals of the Hardy-Stokes functions and it is found that the relevant 
dynamical field parameter is = ^v"^- The authors find that in the region of parameter 
values where x is the smallest parameter in the problem, ^ or eBp± <^ niim'^, 
the self-energy has a dependence on pj_ only through a negligibly small imaginary part 
proportional to e^^^^^. A calculation of the contribution S" to the neutrino self-energy 
from the nth charged lepton Landau level [61] (in conjunction with the exact W propagator) 
also shows that, for low neutrino transverse momentum, E" does not depend on p\, since 

which, for <^ m^r, equals one. Therefore 



pj^ only appears in the factor In ^ 1 + ^, 



in the low transverse momentum limit p]_ <^ fn^, the self-energy does not depend on 
the parameter ^ and one can neglect it from Eq. fl27|) . Since By/ = m^/e ^ 10^'' G and 
mn <^ mw, one can always take 77 -C 1 and <^ 1. After an integration by parts of the 
factor in front of I obtain 
,2 



+ (A coth z sinh zu — Au + u sinh zu) L + 



u coth z cosh zu — A coth z sinh zu + Au) + 
sinh(z + zu 



u- 



sinh z 



where ^± = ^ — ^11 is used and 



(c.t. 



A 



9' 



'1 — u)Xe, + u 
V 

r°°dz /■! 



+ (c.t.)H/(29) 



(30) 



du e-^^n^L. 

JO Z J Q 



(31) 



I follow the same procedure to manipulate the expression of the bubble diagram with the 
scalar, and obtain the following 

-.2 rl 



{u coth z cosh zu — A coth z sinh zu + Au) ^ {XeL + e^R) + 



cosh(z — zu) 

+ (Acoth2;sinh2;M — An — Msinh2;M)^|| (A^L + e^R) — Xirriu — h 

smhiz — zu) sinh(2; — , r t-,\ 

+Aim^ — 0-3 - u — 0-3^11 [XiL + tyR) 



sinh z 
+ (c.t.)$ 



sinh z sinh z 

where I introduce the parameter e^, = ml/m'^r and the counter-term is given by 

9' 



(32) 



c.t.)$ = TTTT^ / — du e~^^[u^ {XiL + t^R) - A^m^]. 



(33) 



The self-energy operator I obtain, S(p) = Siy(p) + S4.(p), is valid for p_L <^ mvy and has 
the following Lorentz structure 

S(p) = aL^+ + cl(p(^7) L+ aii:/+ + Cij(]?(^7) i? + m;. [Ki + ii^r2(7V^7)] (34) 



smce 



CJ3^||L = -(p^7)L, a3^||i? = {p^^)R 



(35) 



and (73 = — 1(7(/}7). This Lorentz structure is in complete agreement with the findings of 



Refs. m 



17|. 



At this stage, we should discuss the meaning of the coefficients appearing in Eq. (IMI) : a^, 
hi and contain contributions from both diagrams, but the contribution from the diagram 



with the scalar is suppressed by a factor of A^. The coefficients a^, and Ki are completely 
absorbed by the neutrino wave-function and mass renormalization. The b and c coefficients 
are relevant for neutrino dispersion but, to lowest order, the dispersion relation depends 
only on bi jo], thus the most relevant of these coefficients. The two coefficients b^ and cr 
are suppressed by a factor of relative to 6^ and cl, and therefore play a less important 
role. However, in the case of a non-diagonal neutrino mass matrix, b^i and cr might produce 
a modification of neutrino mixing in the presence of a magnetic field. The K2 coefficient, 
along with cl and c/j, is needed for the calculation of the neutrino magnetic moment 17|. 

In the past, several authors have attempted to calculate bi jsl, 0, Q, Isl, [q], but all we know 
so far are its values in the limiting case of a "weak field" eB <^ mj, and a "moderate field" 
TiT'iy- Kuznetsov et al. [Gj found that 6^ — 2471^ 

■T]^ (in Xi - fj in the case of 

weak field and b^ = 1^:;^^'^ (Inr/ — 2.542) for a moderate field. Their results are obtained 
using an expansion for the VT-propagator where only terms up to second order in powers 
of the expansion parameter eB are retained. In this paper 6^ and all the other coefficients 
appearing in Eq. are calculated using the exact ly-propagator, and simple analytic 

forms for these coefficients will be obtained that are very accurate for < eB <^ m^. It 
will be shown that, in the case of weak and moderate fields, the expression obtained here 
for bi agrees with the values obtained in Ref. [Q] in the appropriate limits. 

Eqs. fl34|) and fl29l) indicate that, to calculate 6l, we need to evaluate the following double 
integral 

/oo rl 
dz / du e~^^(A coth z sinh zu — Am + u sinh zu). (36) 
Jo 



9' 



167r2 

We start by evaluating the last term 

/oo rl 
dz / du e~^^u sinh zu. (37) 
Jo 

After an elementary but tedious integration we evaluate J{r], Xg) exactly and find 

= (1 - X,y ^ ,f- + 2(l + ,-A,)^ '° [—) - 2(1 - , - A,)'^ '° [—) ■ 
Since 1] <^ 1 and ^ 1, we can expand in the two small parameters to obtain 

J{r], A,) = r/2[i + 0(A,ln A,)] + 0{r]'). (39) 

To integrate the remaining terms of Eq. (l36l) I do the w-integration first and obtain 

/oo rl 
dz / du e"^'^(A coth z sinh zu - Am) = 'n'^G{r]) + 7]^F{Xe/ri) + rfH{'q, Xi) - r]I{'q) (40) 
Jo 
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where I took 1 — A/> ~ 1 and 2 — A/> ~ 2 and introduced the four functions 



F{x) 







zcothz - 1) 2- 



— (2; coth z — 1) 
Z'^ 



G{x) = 

/CO rj ^ 
^coth^e-^^/" 
Z'^ 



+ xe 



-z/x 



-(2 + -) 

X X 



^1-7] 



2\2 



-2 + 



A z 



Tjl — rf 



A 

z 



and 



J(x) = r'^cothze-'^' 
J z^ 



sinh z + sinh z + 2x cosh z 2; sinh z z cosh 2; 



(l-a;2)2 

A numerical evaluation of F{x) shows that 



+ 



+ 



X — x" 



(41) 

(42) 
(43) 



2(x + z) 

X 



F{x) ~ — ln(l + 6x) 
3 



(44) 



(45) 



with high accuracy. G{x) is evaluated analytically for small x by introducing a regulator 
and using the following series expansion of the hyperbolic cotangent 



1 ^ 1 
cothz = - + 2z}^ - — 

z [n-Ky + z'^ 



to obtain 



Gix) = lim 



Vc(2-e)i?f-,l ' 



'-xT(e] 



- 1. 



(46) 



(47) 



,2 2 

Here C(2 — e) is the Riemann zeta function, i? 1 — |^ is the Euler beta function and r(e) 
is the Euler gamma function. After taking the limit we find 



G{x) = IHtt/x) + l-fE - ^C'(2) - 1 



(48) 



3 ' ' ' 3 TT^ 

where je = 0.5772 is the Euler-Mascheroni constant and C'(2) = —0.9375 is the first deriva- 
tive of the Riemann zeta function. Last we evaluate H{r],Xg) and I{f]). After tedious 
integrations we find 

v'H{r^,Xe) = 0{v'), (49) 

and 

2 

r^I(r^) = -!L + 0{v'). (50) 
Once we insert the expressions of F,G,H,I and J into Eq. (136!) we obtain the following 



9' 2 



167r 



■ln(7r/r7) + ^7i^- AC'(2) 



lnh+6^j + 



6 



(51) 



9 



In the case of a moderate field Xe/r] 1 and F^Xi/rf) = 0, and therefore Eq. (l5Ti) gives the 
following 



247r2^ 



lnr/ + -^C'(2)-ln7r-7s-^ 

TT 4. 



(52) 



where ^C'(2) - logvr - - - 
interesting to notice that — In 6 — 



—2.5418 and it confirms the result of Ref. ^]. It is very 



-2.5418 = ;^C'(2)- log vr- 7^ 



|, and this allows 



us to simplify significantly Eq. ( l5Ti) and write it as 



247r2^ 



In ( ^ + A, 



(53) 



This analytic expression is valid for ^ cB <^ and agrees, as I have already shown, 
with the known value of Bl in the moderate field limit. It also agrees with the weak field 
value of bi obtained in Ref. jo], since in the weak field limit rj <^ and therefore = 
{in — I j . An exact numerical computation of 6l has also been done, and its results 



247r2 '/ 



are reported in Figure 1, where the exact value of x is shown for < < 1 and for 
the three neutrino species. Notice that hi diverges as rj approaches one. 

All other coefficients of Eq. ( l34l) have been evaluated and are listed below to leading 
order 

ai -- 



9' 2 



487r' 



In h§r, + A, 



bR 



487r' 



cl 



ln^- + A. 



+ 



327r2 

2 



Cr 
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647r2 



^1 = , 



647r2 



967r2 



A^r7(lnA^ + l) 



(54) 
(55) 
(56) 
(57) 
(58) 
(59) 
(60) 



The following 



i coth , 



1 



~zX/ri 



X 

z—e 

V 



zX/ri 



(61) 



z^ \ zj\ 7] I o \ \/eri ) 

las been used when evaluating a^. The expression for cl agrees with that obtained in Ref. 
gI, the expressions for the other coefficients have not appeared in the literature. 
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IV. DISCUSSION AND CONCLUSIONS 



I have calculated the neutrino self-energy in a magnetic field to one-loop order using 
the minimal extension of the Standard Model with massive Dirac neutrinos. My results 
for all the invariant coefficients of the self-energy are valid for < eB ^ and for 
neutrino energies < E <^ niw and are reported in Eqs. fl53H60l) . These results make the 
distinction between weak and moderate fields obsolete by providing simple and very accurate 
analytic expressions for all the coefficients, and allow us to evaluate the neutrino self-energy 
around the critical magnetic field value = m^/e ~ 4.41 x 10^^ G, where none of the 
previous calculations are valid. Some of the implications of these findings are discussed in 
the remaining part of this section. 

The Dirac equation for left-handed massless neutrino is 

[^-S(p)]Lz/,(p) = 0, (62) 

and once we insert into it the expression of the self-energy with rrii, = 0, we find the 
inverse neutrino propagator 

= ^ - ai^ -hL{p(p(p'^) - cl{p^'^). (63) 

By squaring the inverse propagator and setting it equal to zero we obtain the dispersion 
relation for massless neutrinos 

(1 - ai)y - {2hL -bl-cl- 2hLaL) p] = (64) 

which implies 

p = l+ (^feL + aL6L-y-yjsinV (65) 

where E and p are the neutrino energy and momentum and is the angle between B and 
p. Since ctL&L, h\ and c\ are of higher order, in a perturbative sense, than they can be 
neglected and, using the coefficient hi obtained in this paper (l53l) . the dispersion relation 
can be written as 



IpI ~ 372772 m 



2 

w 



In 



\ 3 



eB I ^2 



+ mi 4 



sin 0, (66) 



where the term with the bracket is the magnetic field contribution. Eq. (166|) is valid for 
neutrino energies E <^ mw and applies to both z/^ and z/^. 
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The neutrino dispersion relation in a magnetized medium has been studied in previous 
papers . In particular, a CP-symmetric medium with temperature rUe <^ T <^ mw 

and magnetic field eB < was studied, because these conditions represent reasonably 
well the early universe plasma between the QCD phase transition and nucleosynthesis. The 
following dispersion relation 



IpI 3 



' ' ' ' cos 6 + -—: — ^ sin 6 m 



15 \Tn\ "^w J "^w STT^m^ ml 



(67) 



was derived for z/g and Ve-, where the first term is the pure plasma contribution [20| and 
the other two terms are caused by the combined infiuence of plasma and magnetic field. 
Eq. (1671) was obtained under the assumption that the magnetic field induced pure vacuum 
modification of the neutrino dispersion relation was negligible. It turns out that the vacuum 
modification could be as large as the second or third term in Eq. fl67|) and, once we include 
the vacuum modification presented in Eq. (!66|) of this paper, the dispersion relation for z/g 
and Vf. in a CP-symmetric plasma with magnetic field becomes 

l-K^T^ (I 2 \ T^eB 

~ — — cos (f)+ 

15 Vm| mfyj 



E 

1 + 



IpI 3 



^ ■ sin^ 6 f In — - — In 



(68) 



and is valid for rrie ^ T ^ mw and < eB < T^. When eB ~ the pure vacuum 
modification can be as large as other terms of Eq. (l68l) and must be included in the neu- 
trino dispersion relation. This dispersion relation leads to an anisotropic neutrino index of 
refraction, causing neutrinos that move in the direction of the field to feel more the effect of 
the magnetized plasma (second term inside the bracket), while neutrinos moving perpendic- 
ularly to the field feel more the magnetization of the vacuum, similarly to what was found 
for strong fields in j^. 

Another interesting application of my result, for the case of massive neutrinos, is the 
resonance enhancement of neutrino oscillations of the type i/g ^ u^t. While it is well known 
that the medium alone can produce neutrino oscillation (i.e. the MSW effect), I want to 
explore here the role of magnetic fields in neutrino oscillation enhancement. Conditions 
for resonance enhancement could be present, for example, inside an exploding supernova 
if a strong magnetic field is generated inside the exploding star, allowing for more energy 
transferred to the stellar matter by the u^. The mixing angle 6b in a magnetized medium is 
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determined by the following relation 

sin^ 26 B = o (69) 

\Amf,\ J 

where Ami ^ squared- mass splitting and vacuum mixing angle in the z/g, i/i 

system and Ve and Vi are the z/g and Ui effective potentials in the magnetized medium. The 
plus sign applies for rrii < mi and the minus sign for m^ > mi, and z = 2,3 for £ = /i, r. 
Even if the vacuum mixing angle 6 is very small, the mixing angle in magnetized matter is 
9b = ^ ii mi > mi and — is positive and satisfies the resonance condition 

V,-Ve=^-^^cos29. (70) 

First I consider neutrinos propagating through a magnetized charged medium, where the 
chemical potential that displays the asymmetry between particles and antiparticles is ;U 7^ 0. 
In the case of a weak magnetic field eB ^ yU^, the pure magnetic field contribution to the 
effective potential is obtained immediately from Eq. ( l66l) and, once we include it into Ve — Ve, 
Eq. fITOl) becomes 

where N^. is the electron number density and T is the temperature of the medium. The first 
term of Eq. (17T]1 is the thermal contribution of the magnetized medium to Ve — Ve and was 
obtained in Ref.Q] for a nonrelativistic and nondegenerate electron gas. 

For a strong field eB ^ fi^ the thermal contribution of the magnetized medium to the 
neutrino self-energy was calculated in Ref. [Oj and, once we include it into Ve — Ve, the 
resonance condition becomes 

y2G.e-i/--(iV0 - iV°) (1 + cos 0) - -^^^E sin^ In ( '4±^] = ^ cos 26 

(72) 

where and Ng are the electron and positron number densities in the lowest Landau level. 
Using a supernova core density of 10^^ Kg/m^ we have \/2GpNe — 3.8 eV and, neglecting 
the neutrino masses and for = 10 MeV, I find that a magnetic field strength B ~ 10^^ G 
is required for the resonance transition Ve ^ to occur inside an exploding supernova, far 
exceeding the magnetic field believed to exist inside the supernova. 

The case of a magnetized neutral medium (/^ = 0) should also be explored. It was shown 
in Ref. 2l| that, in the case of a strong magnetic field eB ^ T^, pure magnetic neutrino 
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oscillations are possible in a CP-symmetric magnetized medium. I find that, for ^ eB 



and ^ m^, the resonance condition for pure magnetic 



oscillations is 



Gf {eBf 



■E sin^ 



TT^ m 



1 

m,, 
ln-| 

mi 



In 



eB + 6m^ 



2 \ n 



|Am|] 



cos 20 



12- 



(73) 



Since Am^i = 8x 10"^ eV^ and cos 2^12 = 0.51 [22], Eq. fl75]) shows that for a neutrino energy 
E = 1 GeV resonance will occur when B ~ 2ml/ e — 8-8 x 10^^ G. This finding can have 
interesting implications for cosmology, showing that resonant neutrino flavor oscillations 
could be caused by a moderate magnetic field in the primeval plasma. 

The magnetic field induced resonance transition i/^ i/^ in a magnetized charged medium 
should also be investigated in this context and, for this transition to occur, the resonance 
condition is 



2E 



{E sin^ 0) In 



eB + Qmi 



eB + 6m?, 



0, 



(74) 



3V27r2 m^ 

since the charged current contribution to the effective potential of and in a medium 
is absent. Eq. (17^ shows that resonance could occur only if Am|2_is negative, which is not 
ruled out [22]. For E = 10 MeV and Am^g - -2-5 x 10"^ eV^ [22] I find that a magnetic 
field B ~ 8.5 x 10^ x v^cos2^i?e is necessary for resonance to occur. Since the accepted 
value of the mixing angle is 37° < < 45° 2^ and B^. ~ 4.41 x 10^^ G, we will have 



resonance conditions for B < 2 x 10^^ G, a magnetic field strength that could exist inside 
magnet ars. 

One more case where the results obtained in this paper are relevant, is the calculation 
of the magnetic moment fi,y^ of Ui within the minimally extended standard model of the 
electroweak interactions containing an S't/(2)-singlet right-handed neutrino. The neutrino 
magnetic moment can be written in terms of the self-energy coefficients of Eq. fl34|) as 17| 



m,. 



where 



cl 



9' 



167r2 

K2 



dz / du 



sinh z 



-u 



sinh(z + zu) — — sinh(2; — zu) 



167r2 



A, 



dz / du e 



-zA 



sinh(z — zu) 
sinh z 



(75) 

(76) 
(77) 



14 



and cr is given by Eq. (l58ll . Once we retain all the sub-leading A^-corrections to cl and K2, 
we have 

g'^r] 1 



17 



A 



3 A^ + 5A; + — + 2A^ In \ - 5A| In A J + 0{n) (78) 



327r2(l-A^)3 V' 2 ■ " ■ 2 
g'^r] 1 /A 



327r2 (I - XfY \ 2 



(lnA, + l-A,) + 0(r/^) 



and, for m^^ <^ m^, we find 



(0) 



7 1 
(1 - A,)3 " 2^' + " ~ 2^' 



where the leading term is [23|, |24 1 



_ 3eGFm^ 



(79) 



(80) 



(81) 



The neutrino magnetic moment of Eq. (IHOll agrees with the results of Refs. 25|, l26|], obtained 
by different methods. 



Acknowledgments 

A. Erdas wishes to thank Marcello Lissia for helpful discussions, the Department of 
Physics of the University of Cagliari and the I. N.F.N. Sezione di Cagliari for their contin- 
ued support and the High Energy Theory Group of the Johns Hopkins University for the 
hospitality extended to him during his several visits. 



[1] G. G. Raffelt, Stars as Laboratories for Fundamental Physics (University of Chicago Press, 
Chicago, 1996). 

[2] A. Erdas, C. W. Kim, and T. H. Lee, Phys. Rev. D 58, 085016 (1998). 
[3] J. C. D'OHvo, J. F. Nieves and P. B. Pal, Phys. Rev. D 40, 3679 (1989). 
[4] P. Elmfors, D. Grasso and G. Raffelt, Nucl. Phys. B 479, 3 (1996). 
[5] A. Erdas and G. Feldman, Nucl. Phys. B 343, 597 (1990). 

[6] A. V. Kuznetsov, N. V. Mikheev, G. G. Raffelt and L. A. Vassilevskaya, Phys. Rev. D 73, 

023001 (2006). 
[7] G. McKeon, Phys. Rev. D 24, 2744 (1981). 

15 



[8] E. Elizalde, E. J. Ferrer and V. de la Incera, Annals Phys. 295, 33 (2002). 

[9] E. Elizalde, E. J. Ferrer, and V. de la Incera Phys. Rev. D 70, 043012 (2004). 

[10] R. C. Duncan and C. Thompson, Astrophys. J. 392, L9 (1992). 

[11] C. Thompson and R. C. Duncan, Mon. Not. Roy. Astron. Soc. 275, 255 (1995). 

[12] C. Thompson and R. C. Duncan, Astrophys. J. 473, 322 (1996). 

[13] A. Brandenburg, K. Enqvist and P. Olesen, Phys. Rev. D 54, 1291 (1996). 

[14] M. Joyce and M. E. Shaposhnikov, Phys. Rev. Lett. 79, 1193 (1997). 

[15] D. Grasso and H. R. Rubinstein, Phys. Rept. 348, 163 (2001). 

[16] J. Schwinger, Phys. Rev. 82, 664 (1951). 



[17] A. V. Kuznetsov and N. V. Mikheev, hep-ph/0605114 



[18] A. V. Kuznetsov and N. V. Mikheev, Phys. Atom. Nucl. 70, 1258 (2007). 

[19] W. Dittrich and M. Renter, Effective Lagrangians in Quantum Electrodynamics, Lecture Notes 

in Physics Vol. 220 (Springer- Verlag, Berlin, 1985). 
[20] D. Notzold and G. Raffelt, Nucl. Phys. B 307, 924 (1988). 
[21] E. J. Ferrer and V. de la Incera, Int. J. Mod. Phys. A 19, 5385 (2004). 
[22] C. Amsler et al. (Particle Data Group), Phys. Lett. B667, 1 (2008). 
[23] B. W. Lee and R. E. Shrock, Phys. Rev. D 16, 1444 (1977). 
[24] K. Fujikawa and R. E. Shrock, Phys. Rev. Lett. 45, 963 (1980). 

[25] L. G. Cabral-Rosetti, J. Bernabeu, J. Vidal and A. Zepeda, Eur. Phys. J. C 12, 633 (2000). 
[26] M. Dvornikov and A. Studenikin, Phys. Rev. D 69, 073001 (2004). 



16 








-14 -12 -10 -8 -6 -4 -2 

Log 7/ 

FIG. 1: Exact value of br x for 10^^^ < < 10*^ for the three neutrino species. For 
—oo< logry < — 1 the analytic expression of 6^ obtained in Eq. ([53l) of this paper is in within 
0.27 % or less of the exact numerical value. 
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